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CONVERGENCE ACCELERATION PROCEDURE 

FOR THE 

METHOD OF STEEPEST DESCENT 


ABSTRACT 


A procedure is proposed which accelerates the conver- 
gence rate of the steepest descent or gradient optimization 
methods. The previously suggested procedures of selecting 
a preferred gradient step size for each iteration is extended 
by defining an easily determined, time dependent weighting 
matrix that approximately extremizes the penalty function or 
performance index. Numerical results with this modification 
are obtained and compared with results obtained by applying 
the conventional technique. A significant acceleration in 
the shaping of the optimal control program is realized. 

INTRODUCTION 

In studying procedures for accelerating the convergence 
rates of the classical gradient methods, it is necessary to 
first define the optimization problem. The purpose of this 
investigation is stated, and a brief background sketch is 
made . 


Definition of the Optimization Problem 

In one class of optimization problems, in particular 
the spacecraft trajectory optimization problem, it is desired 
to determine the history of the control variables in such a 
manner that certain specified initial and terminal constraints 
are satisfied while some performance index is extremized. 

The control variables are unspecified inputs to the system 
which may be chosen to control the spacecraft state; i.e., 
the position and velocity. The initial and terminal constraints 
are simply conditions on the positions and velocities that 
must be satisfied at the initial and terminal time, respectively. 
The performance index is usually a scalar function associated 
with the spacecraft performance and is the quantity to be 
extremized. 

The terminal constraints are handled in either the so- 
called ’’hard" or "soft" form. In the "hard" form an effort 
is made to satisfy tlu terminal constraints identically while 
in the "soft" form the constraints are satisfied only approx- 
imately. It is in the former case that the performance index 
approach is taken because this index is extremized separate 
from the satisfaction of the terminal constraints. It is 
with the latter case that the penalty function concept 
emanates; i.e., a certain penalty is accepted because of the 
less stringent demand of only approximate terminal constraint 
satisfaction. 

Purpose of the Investigation 

The ultimate purpose of this investigation is to develop 
an insight into the convergence characteristics of some of the 
direct optimization methods. This ultimate purpose is approached 
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by satisfying the following secondary objectives: 

(1) Increase the understanding of the currently 
popular optimization methods so that the deficient 
convergence characteristics of each method are 
discovered . 

( 2 ) Extend and modify these methods to eliminate the 
deficiencies . 

(3) Formulate and successfully implement a realistic 
example . 

( A ) Compare the convergence characteristics of the 
proposed procedures with those derived from 
previously proposed schemes. 

Background Study of the Gradient Methods 

An analytical development of a trajectory optimization 
theory was published by Kelley^ ^ in I960. This method, 
referred to as the gradient method, is based on an extension 
of some ideas presented by Courant in 19^1. A similar formu- 
lation was made, simultaneously and independently, by Bryson, 
Denham, Carroll and Mikami^ 2 ^, and Bryson and Denhan/ 3 \ 

( 4 ) 

Kelley, Kopp, and Moyer presented an analysis of several 
gradient techniques using inequality constraints on the 
control variables and a penalty function concept for handling 
terminal constraints. In an effort to determine the thrust 
steering program for the optimization of a second stage 
booster, Pfeiffer^ 5 ^ developed a method of "critical direction" 
which is similar to the gradient techniques of Kelley and 
Bryson. 
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In 1963» more attention began to center around conver- 
gence acceleration for the gradient methods. Wagner and 
Jazwinski^ 6 ^ presented a gradient method incorporating both 
terminal and instantaneous inequality constraints. This 
investigation also included an interesting method for deter- 
mining the control step size magnitude that should be taken 
in the gradient direction to approximately maximize the 
decrease in the penalty function. A new step size is calcu- 
lated for each iteration. This scheme involves making three 
trial forward integrations with different control step sizes, 
and recording the three resulting penalty function values. 

A second order polynominal is fitted through these points, 
and the step size that corresponds to the minimum value of 
the penalty function is selected for the next iteration. This 
method, therefore, takes full advantage of each adjoint inte- 
gration by selecting an optimal step size for that particular 
iteration. 

Rosenbaum , also in 1963, developed a method similar 
to a closed-loop guidance scheme that "provides rapid con- 
vergence for a variety of missions." The distinctive feature 
of this method is that the control step size in the gradient 
direction is calculated and becomes a time dependent quantity. 
The significant result is that larger deviations from the 
nominal trajectory can be tolerated while still satisfying 
the terminal constraints; thus, it is possible to move more 
rapidly toward the optimal trajectory. The approach is 

similar to the A-matrix control scheme proposed by Bryson 
( 8 ) 

and Denham' . Unfortunately, the rates of convergence relative 

ft 

to previously proposed methods are not adequately illustrated. 

In 1964, Stancil proposed a slightly different approach 
to the inherent gradient convergence problem. This approach 
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( 7 ) 

is similar to Rosenbaum in that a time dependent weighting 

matrix is calculated. Basically, the formulation follows a 

suggestion made, but not used, by Bryson, Denham, Carroll and 
( 2 ) 

Mlkami , in which the current control program is averaged 
with the Kulerlan control. The procedure eliminates the guessing 
of the performance index decrease and other weighting factors. 
Again, the convergence characteristics of this method are 
not directly compared with the previously proposed techniques. 

Lewallen^ 1 0 ^ , in 1966, derives a time dependent and 
easily determined weighting matrix which is applicable to 
either minimizing a performance index or a penalty function. 

An analysis and comparison is made using both the proposed 
weighting matrix and the unity matrix. When the proposed 
matrix is used a significant convergence acceleration is 
realized . 


FORMULATION 

The gradient method is formulated with both hard and 
soft constraints even though the formulation, perhaps in 
different form, is presented in available literature. Inclu- 
sion of this information in the present report is encouraged 
(1) to make the report more self contained and (2) to provide 
a basis for the extensions required for convergence accelera- 
tion discussions. 

Gradient Method Using Hard Constraints 
It is desired to determine the control program u(t), 

ft 

where u is an m vector, which will yield an extreme 
value of some performance index 

<f> = $(x f ,t f ) (1) 
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subject to the differential equations of motion 


x *= f(x,u,t) ( 2 ) 

where x is an n vector and u is an m vector, while 
satisfying the terminal constraint relations 

♦ * <l'(x f ,t f ) ■ 0 (3) 

in hard form, where is a q vector. One of the terminal 
constraint relations may be selected as a stopping condition 
for the integration process, 

P = fl(x f ,t f ) = 0 . * (*0 

If the differential equations (2) are linearized about 
some nominal path, the resulting equations become, 

6 x = f x 6 x + f u 6 u ( 5 ) 

where f and f are partial derivatives of f with 

respect to x and u , respectively, and are evaluated on 
the nominal trajectory. 

The equations adjoint to (5) are 

x « -tlx (6) 

where A is an n vector of adjoint variables. This 
equation may be combined with ( 5 ) to yield 

jr-(X T «x) - X T f u «U . (7) 
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Integrating this equation yields 


(A T 6x) f *= J A T f u 6u dt + (A T 6 x) q (8) 

t o 

which is designated the Fundamental Guidance Equation. The 
object now is to determine how initial state variations and 
integrated control variations influence the performance 
index, stopping condition, and the terminal constraints. If, 
on separate trials, the terminal values of the adjoint vari- 
ables are set equal to 




riii x T (t ) = m 

UxL V f' ldx\ 


A ft ( t f ) 


Taft' 

Lax. 


( 9 ) 


where A is an n vector, A is a n x q matrix and 


♦ 




A fi is an n vector, the desired relations are seen to be 


d 


f f T 

- J f U 


T 


* 6u dt + (A 6x) + $ dt r 

<p U 4' O I 


( 10 ) 


dij> 


A* f 

J xj f u 6u dt + Uj6x) 0 + i dt f 

t _ 


( 11 ) 


an 


/ 


x£ f u 6u dt + U„«x) 0 + a dt f 


( 12 ) 
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where 


and d() 


[6() + 0) dt] 



f * 


This formulation allows the specification of an allowable 
step size to be taken in control space defined by 


dS 



1 t 

2 6u W 6u dt 


(13) 


where the step is a weighted quadratic function of the con- 
trol deviation. The weighting matrix W is included to 
improve the convergence characteristics by giving more weight 
to regions of low sensitivity. However, it is often chosen 
unity because of the lack of knowledge concerning the region 
of sensitivity. The criteria used for determining the best 
elements for this weighting matrix are not easy to determine 
and are usually found through trial and error procedures. 

The stopping condition 00 is to be identically satisfied 
so d ft in (12) is equated to zero. The terminal time varia- 
tion dtj. is eliminated from (10) and (11) to yield 


d 4> 


-/ 


X Ifi f u 6u dt + (X Jn 5x >o 


( 1 * 1 ) 


d<p* 


■ r 


X I« f u 5u dfc + (x Ja* x) 


(15) 


where X 




A. - - A 
♦ fl 0 


X , = A . - x_ - 

Ip £ 
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The total variation of the performance index may be 
represented by 


d<$> 


/ 


f 6u dt + (X^fix) 
u o 


+ v 


d^ 


< 


X^’ f 6u dt 


<*Jn {x) o 


+ u 


dS 


- 


j 6u^ V/ 6u dt 


( 16 ) 


where the terminal constraints and the control step are ad- 

T 

joined by the use of the v and p Lagrange multipliers, 
respect ive] y . Since it is desired to determine the control 
variation which corresponds to the maximum change in the 
performance index, the first variation of (]6) must vanish; 
therefore 


/ f 

(x In f u - vT V f u " “ 6uT w) 4 u dt h o . 

^O 

( 17 ) 

This implies that the desired control variation is 


» 


6u = 



u 


(X 


<f>n 



( 18 ) 
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and when this equation is substituted back into (13) and 
(15) » the values of v and p are seen to be 


v 




and 


T - 1 

I - I I I 

<p< j) M M M 

dS - dB T I”' dB 


where 


dB = d* - (X i(in «x) 0 






T 


.-1 


x; n r it w f* x . dt 

ipn u u ijjft 


**♦ - V x l* f u W_1 r u V dfc 

l o 

- r »j c f u w_i f u So at 

l o 


and I ^ is a q*q matrix, I is a q vector, and 
is a scalar. 


Now combining (18) through (23) yields the desired 
control program 


( 19 ) 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 

I 4'4> 


10 


6u 



f T (* 

U V 


X I ” 1 I ) 

<pn ipip lp<r 


dS - dg T I"’ de 

»» 

T 

I - I I I 

<p <p ip 4> ip ip ip <p J 


+ W 


“ 1 x I” 1 

u i i>n (pip 


d 3 


( 2*0 


where the positive (negative) sign is used if <p is to he 
maximized (minimized). The previous control program is now 
modified by 


u 

new 



6u . 


(25) 


The computational procedure for the Gradient Method 
using hard constraints may be summarized as follows: 

(1) Integrate the n differential equations of motion 

(2) forward, using an assumed control program and 
the desired initial conditions. This integration 
is continued until the stopping condition (*0 is 
satisfied. The state variable values are stored 
at each integration step. 

(2) Integrate the adjoint equations (6) backwards 
q+2 times with the starting conditions (9)* The 
coefficient matrix f is formed from the state 

A 

variables stored during the forward integration. 

(3) Integrate the I equations (21) through (23) back- 
wards simultaneously with the adjoint equations 
using initial conditions of zero to yield values at 

t for I.., I, . , and I,. . 
o » ipip ip<p * 

(4) Select a desired improvement in the terminal dissatis- 
faction dip for the next iteration. 


(5) Select a reasonable value for the mean square 
control deviation from the previous control 
program by using 


dS 


\ 

2 


6u 2 

ave 





which provides a value of dS . 

(6) Use the selected value of dij/ and dS to calcu- 

late the numerator under the radical in (2*0. If 
this quantity is negative, determine the that 

makes the quantity vanish. If it is positive, use 
the quantity as it is. 

(7) Calculate the 6 u as given by (2*0 and alter the 
assumed control program. The quantity dS must 
be decreased according to some selected criteria 
to prevent stepping across the optimal point into 
a nonoptimal region. 

(6) The procedure is continued until the control varia 
tions are less than some preselected value. 
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Convergence Acceleration Using Hard Constraints 

A primary objective of the present investigation is 
to develop an iterative scheme that reduces sonic of the 
arbitrariness and increases the convergence rate of the 
Gradient Method when using hard constraints. Since the 
weighting matrix W , Introduced in (13) is arbitrary, 
some rational basis for its selection is needed. 

The problem is approached by examining an integral form 
of the Weierstrass E-Function which approximates the change 
in the performance index. This change is 


d<$> 



( 26 ) 


where E is the Weierstrass E-Function as developed by 
Gelfand and Fomin 1 1 ^ . The E-Function is defined as 


E 


ft 

f(x , X , t ) - f ( X , X ,t) - 


“jr( X , X ,t) (X - X ) 

ax 


(27) 


and for the system to be considered 

f(x,x,t) *= H(x,u,t) - X^x . (28) 

The asterisks refer to the optimal path, and the absence of 
asterisks refer to any nearby path. From the calculus of 
variations, a necessary condition for the existence of a 
minimum valued performance index is that E be non-negative 
during the interval t Q < t < t f . 
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It is noted, by examining (?) , that a variation in con- 
trol is accompanied by a variation in x , and that a state 
variation will occur only after a finite duration of time. 
Hence, the expansion of (26) is made by considering that the 
control deviation is not accompanied by a change in state. 
The relation (26) is now written 


d<> 



(H - H* ) dt . 


(29) 


The first term In the integrand may be expanded in a Taylor’s 
series about the optimal path at each poinu in time to yield 


^ # # l T * 

H » H + H tl 6 u + 7 6 u 1 H 5u + 

U 2 UU 


(30) 


and substituting the above equation into ( 29 ) and recalling 

* 

that H *=0 on the optimal path results in, 


d<f> 



6u dt . 


(31) 


This equation represents the performance Index change 

associated with the deviation of the control program from an 

* 

optimal control program. It must be stated that H uu is not 
known until the optimal trajectory is converged upon, but the 
expression (31) becomes increasingly accurate as convergence 
progresses. It is during this terminal phase of convergence 
that the Gradient Methods have the greatest need for con- 
vergence acceleration. 


The convergence acceleration technique proposed Jn the 
present investigation uses the expression (31) to approximate 
the performance index change rather than the one previously 
mentioned in (1*0. When this is done and the control 3tep 
size constraint is not included, d$ may be written 


d 4> 



6u dt 


+ 



dt 



6u dt 


(3?) 


Requiring that 6(d$) = 0 yields 


6u 


(H 





v 


and when this equation is substituted into (15) > 


v 



d 6 


(33) 


(3*0 


where 




/' 


qi # i T 1 

X, 0 f (H )“ 1 f A 10 dt 

U UU U 


(35) 


Therefore, the desired control deviation becomes 


su “ (H uu r ’ f u r ;J dB • 


(36) 


15 


In comparing equations (36) and (2*1), it should be noted 

that (36) is simply the last term of (2*1) where W has been 

11 

replaced with H yu . Since the step size constraint was 
eliminated, the control variation can bo controlled by 
requesting only a percentage of the terminal dissatisfaction 
to be corrected by 


dB * -C d0 , 


(37) 


where 0 < C < 1.0 . 

It is interesting to note that the control variation 

law ( 36 ) proposed in the present investigation is similar 

( 1 ? ) 

to the one successfully used by Tapley and Fowler in a 
closed-loop control scheme. 

Gradient Method Using Soft Constraints 

The theoretical development of the Gradient Method using 
soft terminal constraints is similar to that used for hard 
constraints. The primary difference is that the terminal 
constraints are adjoined to the performance index to form a 
penalty function 


P(x f ,t f ) = W o 4 . 2 (x f ,t f ) + V Wj *|(x f ,t f ) ( 38 ) 

\ 

where the W^'s are weighting constants. If these con- 
stants are sufficiently large, minimizing the penalty func- 
' ion is essentially the same as minimizing the performance 
index while driving the terminal constraints to zero. 
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To determine how this penalty function is related to 
initial state variations and the integrated control varia- 
tions, the Fundamental Guidance Equation (8) is used. 
Selecting the starting conditions for the adjoint equations 
(6) to be 


*p(t f ) 


'ajf 

. sx 'l 




’an" 

. 8x Ji 


( 39 ) 


where Xp is an n vector and X fi is a scalar, yields 


dP 


/ 


Xp f 6u dt + ( X p 6 x ) Q + P dtp 


( 40 ) 


do 


/ 


xl f. 6u dt + (x^6x) + n dt ~ . 

U u W O I 


( 41 ) 


If the stopping condition dft is identically satisfied, 
the penalty function change may be expressed as 


dP 




6x ) 


o 


where 



( 42 ) 


( 43 ) 


\ 


17 


Now, it is desired to determine the control variation 
which maximized the penalty function change dP . Adjoining 
an unweighted control step constraint to the penalty function 
change yields 


dP - 


/ f 

X Pf) f u 4u dt 
t 


+ 



6u dt 


+ (*£ fi «x) o . ('I'D 


Requiring that 6(dP) vanish implies that 


6 u 


K f u X P« 


K H 


u 


(^ 5 ) 


is defined as 


where K is a constant equal to 1/p , H 

rji U 

X p Q f . This equation is similar to the one developed by 
Wagner and Jazwinski v . The constant K can be Interpreted 
as a control step size in the gradient direction. 

The penalty function change is evaluated by substituting 
(*I5) into (*I2) to yield 


dP = K 


A T 

/ H H 1 

J u u 


dt 


( 46 ) 


» 
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The computational procedure for the Gradient Method 
using soft constraints may be summarized as follows: 


(1) Integrate the n differential equations of 

motion (2) forward using an assumed control program 
and the desired initial conditions. This inte- 
gration Is continued until the stopping condition 
(*0 is satisfied. The state variable values are 
stored at each integration step. 


( 2 ) 


(3) 


Integrate the n adjoint equations (6) backward 
one time with the starting condition (*I3). The 
coefficient matrix f is formed from the state 
variables stored during the forward integration. 

Having obtained the solution X p ( t ) , the term 

rji rB 

H u = Xp^ f may be formed. The square of H 
may be integrated from t to tp and the step 
size K may be determined by specifying a desired 
penalty function change dP . 


( *1 ) The control variation may be determined from (*I5) 
and applied to the previous control program. 


(5) The procedure is continued until the control vari- 
ations are less than some preselected value. 


It must be noted that the specified penalty function 
change, and hence the step size K , is arbitrary, and 
the judicious selection of K becomes a ke> factor in 
increasing the convergence rate. An automatic procedure for 
its selection is desired. 
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Convergence Acceleration Using Soft Constraints 

The formulation for determining a convergence accel- 
eration procedure when using soft constraints is similar to 
that when hard constraints are used. The performance index 
change or penalty function change in the case of soft termi- 
nal constraints is approximated by 



(** 7 ) 


In the same manner as discussed for hard constraints, this 
relation may be reduced to 


dP ~ 



('18) 


which states that the penalty function change may be 
approximated by a time integral of the Hamiltonian deviation 
from the optimal value. 

Now, the second term in the integrand of (*l8) may be 
expanded in a Taylor’s series about the current (nonoptimal) 
path at each point in time to yield 

H* = H + H 6u + i 6 u T H 6u + . (*»9) 

U a uu 
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Substituting this expression into (*18) results in 


dP 



(H 6u + - 6u T H 6u ) dt 

U 2 UU 


( 50 ) 


which states the penalty function change in terms of the 
control deviations. 

It is desired to determine a control deviation which 
will maximize the penalty function change on each iteration, 
and hence a necessary condition is that 6(dP) = 0 which 
leads to 


6u 


-H 


-l 

uu 



(51) 


where H uu and H are evaluated on the current trajectory. 
This equation implies that the optimal control is in the 
negative gradient direction, weighted by H~ u . By comparing 
(51) with (*15), it is seen that the constant K , calculated 
only once for each iteration, is replaced by the time 
dependent and easily calculated weighting matrix . 


APPLICATION AND RESULTS 


The theoretical developments made in the previous section 
are applied to a realistic example which is difficult enough 
to demonstrate the convergence advantages but simple enough 
to be easily implemented. The example chosen is the two- 
dimensional minimum time, constant low thrust Earth-Mars 
transfer trajectory. The transfer is assumed to leave the 
Earth region with initial conditions corresponding to that 
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of Earth, transfer through heliocentric space to match con- 
dition of the Mars orbit. The orbital parameters used are 
shown in Appendix A, and the differential equations used are 
shown in Appendix B. An example procedure for setting up the 
Method of Steepest Descent using hard constraints Is shown in 
Appendix C. 

The Method of Steepest Descent using hard constraints is 

selected to illustrate the convergence acceleration procedure. 

« 

The value of using the H matrix in the expression for 
control deviation (36) rather than using the unity matrix 
in the conventional formulation (2*1) is illustrated in 
Figures 1, 2, and 3. The three figures represent cases 
where three widely different initially assumed solutions or 
control programs are used. The plots are of thrust angle 
above local horizontal in degrees as a function of mission 
time in days. Each figure also includes the Eulerian or 
optimal control program solution so that the state of con- 
vergence of the other illustrated solutions may be assessed. 
The two remaining solutions on each figure are the ones that 
have been developed after 13 iterations. The curves marked 

by W - I used the conventional technique described by (2*0 

* 

and the ones marked by W = H uu used the proposed technique 
(36). 

The significant fact illustrated is that for all three 
assumed solutions, the control program, after 13 iterations, 
that uses the proposed acceleration procedure is well ahead 
of the conventional procedure in shaping the curve. Both 
procedures ultimately approach the Eulerian program, 
therefore illustration of the converged solution is not 
instructive . 


22 



9-[3uv isruqj, 









C\J 











j°.‘ 

1 tmtrm; 

} C/J 

>* • 






o 



1 







! VO 











•-* • t • • • < 

• 


: © . 

ii 

:.L : _ 

! ° 



3 

* 3 


W t 

— 

II 


:< - 

•v , rnr. 


! (/) 
n 

w V) 
1 U <D 

to c 

^ l 

o o 

o> 

l t< 

P • 

aa>^ w 

JC 

cvj c: 

! rH P 

o 

1 ° 

*0 *H. 

S. <D 

0) P 

i P ?< 

t: cJ 

C ri*H In 

I O 

n <D 

o 

U P 

•> — 

> 

<L> *H 

< 1 ) o 

P 

x: x: 

0) cn 

Eh W 

r a rH 

• • 

0> 

P 

o 

• •»— » rrr — i 





Although it is not shown, for the W *= 1 case illus- 
trated in Figure 2, approximately 12 more iterations are 
required to duplicate the shape obtained by the W * H 
case In only 13 iterations, hence, for this particular 
case, the proposed procedure reduces the computational time 
required to 50 percent of the tine required by the conven- 
tional method. 

One additional piece of information that can be extracted 
from the figures is that o r how the assumed solution influences 
the convergence rate for this particular problem. Since the 
control program state of development is shown for the 13 th 
iteration in each case, a comparison may be made. 

CONCLUSIONS AND RECOMMENDATIONS 

The conclusions of this investigation arc that the 

* 

proposed procedure, where W = H , produces a significant 
acceleration in the convergence rate of the Method of 
Steepest Descent using hard constraints . In one of the cases 
presented, the computational time was reduced to one-half 
of that previously required. A byproduct in the investi- 
gation is seen in that convergence occurs for three widely 
different initially assumed control programs — a real con- 
trast to the highly sensitive indirect methods. It is seen, 
however, that the assumed control program does influence the 
rate of convergence. This is illustrated in the figures by 
comparing the control program shape after 13 iterations for 
each different initially assumed solution. 
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This theoretical development for the Method of Steepest 
Descent with soft constraints has not been verified by 
application to the above example. It 1, recommended that 
this be completed, so that a comparison of the two procedures 
can be made. 
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Appendix A 
ORBITAL PARAMETERS 


Appendix A 
ORBITAL PARAMETERS 


Astronomical Unit, AU 

Orbital Radius of Earth, r 

* e 

Orbital Radius of Mars, v R 

Gravitational Constant 
of Sun, GM 

Initial Spacecraft Mass, m c 
Spacecraft Thrust, T 
Spacecraft Mass Rate, m 


.1^959870 x 10 12 meters 
.10000000 x 10 1 AU 
.15236790 x 10 1 AU 
.1327150*1 x 10 ?1 meters 3 /sec 2 

.67978852 x 10 3 kilograms 
.*10312370 x 10 1 nev;tons 
.10123858 x 10" 14 kilograms/sec 
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Appendix D 

DIFFERENTIAL EQUATIONS 

The differential equations used in this investigation 
(1) the differential equations of motion 


u 


r ? GM + T si n 6 
r “ r 2 m 


Z o = V = - 


uv T cos B 
r m 


z 3 « r = u = f : 


z 4 


- = f 4 

r 4 


(2) the adjoint differential equations 
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EXAMPLE PROCEDURE 


Appendix C 
EXAMPLE PROCEDURE 

The differential equations of motion are integrated 

forward from t with starting conditions 

o 



and some initially assumed control program fl(t). 

The performance index to be minimized is 

* = t f 

and the terminal constraints are 

♦ 1 = u(t f ) » u f = 0 

i>2 = v(t f ) - v f = 0 

^ 3 ■ r(t f ) - r f = 0 . 

The stopping condition used is fl = e(t f ) - o f = 0 , and the 
starting conditions for the backward integration are 







rr 



10 0 0 
0 10 0 
0 0 10 





[0 0 0 1 ]. 


The time rates of change of performance index, 
terminal constraints, and stopping condition are 
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